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In this paper we apply topology optimization to micro- 
structured superhydrophobic surfaces for the first time. 
It has been experimentally observed that a droplet sus- 
pended on a brush of micrometric posts shows a high static 
contact angle and low roll-off angle. To keep the fiuid from 
penetrating the space between the posts, we search for an 
optimal post cross section, that minimizes the vertical dis- 
placement of the liquid-air interface at the base of the drop 
when a pressure difference is applied. Topology optimisa- 
tion proves effective in this framework, showing that posts 
with a branching cross-section are optimal, which is con- 
sistent with several biologic strategies to achieve superhy- 
drophobicity. Through a filtering technique, we can also 
control the characteristic length scale of the optimal de- 
sign, thus obtaining feasible geometries. 
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Superhydrophobicity is a remarkable natural phenomenon, re- 
cently analysed ^ ^ and reproduced artificially ^ ^ by numer- 
ous research groups. Superhydrophobic surfaces show very 
large static contact angles and small roll-off angles for water, 
and these properties are usually associated with self-cleaning 
surfaces. 

A micro- and/or nan o-sca le texture is usually at the origin 
of superhydrophobicityEEl. A drop can reach several differ- 
ent equilibrium states on a textured substrate, as sketched in 
Fig. [T] The effective minimum energy configuration depends 
on the chemical and geometrical properties of the liquid- solid 
interface. We will now focus on superhydrophobicity, which 
is usually associated to the Cassie-Baxter state In this con- 
figuration, the drop is suspended by the protruding features, so 
that its base is in contact with a heterogeneous solid-air sub- 
strate. The apparent static contact angle dcB , according to 
Cassie-Baxter theory, is given by: 



cos dcB = fsl COS 07 - ( 1 - /,/ ) , 



(1) 
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Fig. 1 Sketch of the possible equilibrium positions for a drop on a 
textured surface. The states between orange bracket are accessible 
for hydrophilic materials, those between blue brackets are accessible 
for hydrophobic materials. 



a weighted average between the contact angle for the solid 
substrate (By) and for air (^air = 180^), where fsi represents 
the wetted solid surface per base area of the drop. 
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Fig. 2 A) Interface deformation under applied pressure, for 
hydrophobic materials. If Qtend > ^Y^ the contact line slides along 
the side of the post, in the direction of the yellow arrow B) 
Analogous setup for hydrophilic materials. If Qbend > ^Y, the liquid 
wets the bottom face of the plate, in the direction of the yellow 
arrow. 

If a drop in the Cassie-Baxter state is perturbed, for instance 
if a pressure difference is applied between the drop and the 
environment, the liquid-air interface will bulge, and eventually 
the liquid will begin to flow along the side of the posts when 
the angle dhend (see Fig. [2]) exceeds the contact angle By . 

This effect is particularly important for inherently hy- 
drophilic materials, for which a heterogeneous wetting state 



can be achieved through overhanging structures (see Fig.) 
[2^), even if the global energy minimum will be a Wenzel 
stately (Fig. [T]). Maximising the robustness of the suspended 
drop configuration upon applied pressure is therefore funda- 
mental for effective superhydrophobic surfaces. 

The research about "Cassie mode" superhydrophobicity has 
so far been characterized by a strong dichotomy. On one 
hand, complex hierarchical structures have been fabricated 
and tested experimentally, but their modelling is hard, since 
the structures are usually rough and non periodic. On the 
other hand, there is an active research for the optimal post 
shape to achieve a robust Cassie state, which however usually 
relies on simple shape perturbations to conventional cylindri- 
cal or square posts. In this paper, we take a step in bridging 
this gap, applying the tools of topology optimization. 

Topology optimization^^ is a structural optimization 
method with no intrinsic constraint on the topology of the so- 
lution, which has been applied in such different fields as struc- 
tural mechanics photonic crystal design and microfluidic 
devices We will here apply it to obtain the texture that min- 
imizes the deformation of the liquid-air interface under ap- 
plied pressure, thus making the suspended state as robust as 
possible. We will see that this approach generates interesting 
branching structures, which resemble natural and experimen- 
tally tested superhydrophobic structures. However, the sym- 
metry and length scale of the optimal design can be tuned in 
the numeric optimization procedure, leading to a better under- 
standing and control of such features. 

Modelling and numeric setup 

In this work we will restrict our analysis to a unit cell for a 
square array of posts (Fig|3]), neglecting finite size effects at 
the edge of the drop. We will consider a two dimensional pic- 
ture, in which the liquid-air interface is flat and suspended 
on top of the posts (z=0) in the unperturbed configuration, 
and bulges between the posts to a depth S{x) upon applying 
a pressure difference AP. Such a pressure difference across 
the liquid-air interface can arise for different reasons, such as 
the Laplace pressure due to the drop curvature or the pressure 
upon impact of a drop on the substrate. We also introduce non 
dimensional unit for length /, surface tension a and pressure 
P as follows: 

/ = Lq/, 

Go - 

P = -^P = PqP. 

Here Lq is the characteristic length of the system, which we 
will take as the side of the unit cell (typically few jum), and 



Go can be taken as the surface tension of the liquid considered 
(72.9 mJ/m^ for water at 20 °C). Moreover, since typically 
Lo «lc = yjf^^ where Ic is the capillary length for the liquid 
considered, we can neglect gravity. 

Let us first consider a simple geometry, such as a cylindrical 
post ( cross section is shown in Fig[3p ). The deflection of the 
liquid-air interface among posts can then be described by the 
Young-Laplace equation*^ 




= ^P onD 



< Six) =0 on aOi (3) 

V^(jc)-H = onaD2. 




Fig. 3 A) Sketch of the considered post array. B) Top view of a 
single post cell, for the basic circular cross section. C) Same view of 
a single cell, with a variable cross section (white area) around a 
fixed "nucleus" (red dot) . In the topology optimization procedure, 
the cross section is not fixed but varies according to the field y(x). 
The dashed line in C) shows the reduced computation domain that 
exploits the symmetry of the cell. 

A Dirichlet boundary condition S{x) = is used at the 
boundary of the solid structure 3Di to represent that the in- 
terface is pinned on the ridge of the post. A Von Neumann 
condition VS{x) • H = is applied on the boundary of the unit 
cell 3D2 to account for the symmetry of the post array (in the 
following, we will also exploit the symmetry of the cell to 
work only on one eighth of the domain). 

For the optimization procedure, we will now slightly mod- 
ify this setup. We still consider a "solid" support (Six) = 0, 
red dot in Fig[3p ) in the centre of the domain, but now the 
post cross section is allowed to change around it, in order to 
provide an optimal support for the interface. The distribution 
of material at point x inside the cell is described by the design 
variable y(x), a scalar field which ranges from (completely 
solid) to 1 (completely empty) through intermediate values. 

The field y(x) will be coupled to Eqs|3] leading to the fol- 



lowing formulation of the problem: 

S{x) = 



VS{x)-n = 



AP onD 
on 3Di 
on 3D2. 



(4) 



Where K{y) is defined as: 



{Km 



l)-^-(l-y) 



(5) 



Given the form of Eq. |4]and|5] it is possible to understand 
the effect of the design variable y{x) on the solution S{x). 
Where y{x) = 0, ^(y) is equal to K^ax, which is fixed to be 
a large value. The value AP on the right side of Eq.4 becomes 
then negligible, and the liquid-air interface S{x) does not de- 
form significantly. We therefore recover the "solid" condition 
S{x) 0. On the other hand, If y{x) = 1 ("empty space"), 
^(y) = 1, and we recover the Young-Laplace equation out of 
the support of the post, as in Eqs|3] Intermediate values of 
^(y) do not have a direct physical interpetation, but are re- 
quired for a smooth optimization procedure. The interpolation 
between these two extreme ranges is controlled by the param- 
eter q in Eq. |5] By choosing a sufficiently small value (here 
10~^), it is possible to drive the optimization procedure to give 
a well defined "solid-empty" binary design. 

This formulation also resembles a 2D optimal heat conduc- 
tion problem, where ^(y) corresponds to the distribution of 
conducting material^. 

We eventually need to define an objective function, i.e. a 
quantity whose minimization with respect to y{x) will max- 
imise the support to the liquid-air interface. We choose this 
quantity, called <I> [5'(x),y(x), AP], to be the squared integral 
displacement of the interface (for a given pressure difference 
AP and material distribution y{x)): 

^[S{x),y{x),AP]= [ S^{x)dA. (6) 

JD 

With this choice, we do not control directly the angle be- 
tween the interface and the side of the post, which is indeed 
what would trigger the penetration of the liquid among posts. 
However, Eq. |6]is easy to evaluate through the optimization 
procedure, and its minimization naturally constrains the max- 
imum bending angle of the interface although there might 
be fluctuations along the post ridge. 

At every iteration, the topology optimization code 
changes the value of y{x) over the domain and evaluates 
<I>[5'(x),y(x),AP] and the sensitivity ^^y. We then use this 
information as input to find the configuration of y{x) that min- 
imizes the objective function <I>, using the method of moving 
asymptotes (MMA)^^. Details on the sensitivity analysis and 
the implementation of the code can be found in the paper by 



Olesen et al.^. We will also introduce a constraint on the 
maximum solid fraction per unit cell as: 



\-y(x)dA<fsi. 



(7) 



Remembering Cassie-Baxter relation cos0c5 = /^/ cos By — 
fsi) ^ Eq. [rjcan conveniently be interpreted as a constraint 
on the static contact angle shown by a surface patterned in this 
way. 

The specific coupling ^(y) we use in Eqs. |4]will generate 
a structure connected to the boundary 3Di, which "radiates" 
the support to the y regions This effectively make our 
analysis a shape optimization with many degrees of freedom, 
while the general topology optimization routine we use could 
as easily generate disconnected topologies. 

There are a few reasons for the choice of connected design. 
First, it is well known that dense and thin posts, ideally down 
to the nanometer scale, offer i ncrea singly better support to 
drops in the Cassie-Baxter statel^El. However, it is perhaps 
more interesting to optimize the shape of a single texture el- 
ement, which can then be scaled up or down in size accord- 
ing to fabrication and performance constraints. Second, if we 
are interested in obtaining a hydrophobic behaviour from hy- 
drophilic materials, overhanging structures are required. In 
this perspective, the central support in our optimisation can 
be considered as the stem of the post (see Fig. [2]), while we 
effectively optimise the cross section of the top plate. Even- 
tually, we argue that connected structures would show higher 
mechanical robustness than hair-like features, in particular to 
buckling and shear loads. This latter property is of great rel- 
evance for practical fabrication purposes, since most practi- 
cal application would include significant stresses for the sub- 
strates.^^ 

A final remark regards the length scales in the optimal de- 
sign: at every iteration in the optimisation routine we calculate 
a smoothed version y{x) of the design variable y{x) , applying 
a diffusion stepE^: 



L',,ffy%x)=%x)-y(x) 



(8) 



While calculating the sensitivity, y{x) is then used. This 
process allows to control the minimum size of the features ap- 
pearing in the optimal design. As we will discuss in the next 
section, without filtering small length scale features would ap- 
pear in the optimal design, ideally down to the mesh scale. 
However, these small solid features surrounded by empty 
space are transformed by the diffusion step in a homogeneous 
area with intermediate y{x) value, and thus are penalized by 
the ^(y) function, which favours a binary solid-empty solu- 
tion. The main advantage of this technique is its formulation 
in terms of a partial differential equation, which relies on the 
same numeric tools used for Eqs. |4] 



The actual implementation of our optimization routine uses 
a Matlab code, that relies on the commercial software COM- 
SOL to solve the partial differential equations at every itera- 
tion step. 

Discussion of optimized designs 




Fig. 4 A) Top view of a cylindrical post with solid fraction 
fsi = 0.25. B) Displacement plot for applied pressure AP = Pq- C) 
Optimized material distribution with solid fraction /^/ = 0.25. D) 
Displacement plot for the topology optimized design. 

In the following, the pressure difference acting on the inter- 
face has been fixed as AP = Pq • In Fig. |4]We compare the 
performance of a cylindrical post (A) and an optimized design 
(C) inside a unit cell. The surface plots displayed on the right 
(B-D) show the vertical displacement S{x) obtained through 
Eqs|3] For both structures, the solid fraction is fsi = 0.25. It 
is easy to see the enhanced performance of the topology op- 
timized structure, with the mean displacement reduced by a 
factor 10. It is clear that the branching in the optimal struc- 
ture increases the length of the contact line, where the surface 
tension acts on the side of the post to balance the effect of the 
applied pressure difference. This result in a reduction of the 
interface deformation. However, we think that just choosing 
a meandering cross section would not improve dramatically 
the performance. Lobaton and Salamon^, for instance, con- 
sidered a simpler sinusoidal perturbation to a circular cross 
section. While significantly increasing the contact line length, 
such a shape modification showed modest improvement in the 
critical pressure value. The added feature of our optimal de- 
signs is the convenient placement of the branches, that adjust 
to the cell shape (here a square unit cell, however analogous 
solution have been tested for hexagonal lattices) to reduce the 
size of the gaps between solid features. We therefore argue 
that the significant reduction in the surface displacement arises 
from the interplay of optimal location of the main branches 
and increased contact line length coming from the secondary 
branches. 





Fig. 5 Optimal design for f^i = 0.3 and L^iff =0.5, 1, 3 times the 
meshsize h^esh- The radius of the yellow dot in each column is 
equal to L^iff- 



This physical picture makes it easy to understand the effect 
of the filtering length L^i/f on the optimal design. The designs 
shown in the upper row of Fig. [5] were obtained by solving for 
the domain shown in the bottom row. The yellow dots have a 
radius equal to L^i/f. 
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Fig. 6 Optimal designs as a function of solid fraction /^/. The red 
dot represents the fixed support in the optimisation procedure. In the 
chart, the mean vertical displacement of the liquid-air interface is 
compared for cylindrical posts and topology optimised ones. 

It can be seen that, for any value of L^iff, the structure 
branches along the diagonals of the square cell, thus filling 
the largest gap between two posts. If the resolution is suf- 
ficiently fine, further branching appears, with new branches 
filling the gap among the diagonals. The process continues for 
even smaller length scales and we get an overall quasi-fractal 
behaviour. 

It is possible to see how the filtering procedure constrains 
the minimal length scale in the optimal design. This allows to 



obtain structures suitable for fabrication, i.e. with a feasible 
amount of branching. 

The fractal-like structures resemble several biologic sur- 
faces (such as the lotus leaf), which use analogous (although 
three dimensional) multi- scale structures to achieve their su- 
perhydrophobic properties. A three dimensional optimization 
would be very intensive in terms of computation. It is how- 
ever possible to complement the suggestions from topology 
optimization with general knowledge from superhydrophobic 
surfaces, to get an even more effective texture. Indeed, most 
artificial and natural superhydrophobic surfaces are character- 
ized by a micron scale texture with superimposed nanometric 
roughness. The cross sections shown here should be consid- 
ered an optimal micron scale pattern, over which nano-grass 
features can be grown, thus achieving a multi-layer support 
for the interface (this procedure is currently been considered 
in collaboration with Nis K. Andersen and Rafael Taboryski, 
and will be the subject of a future publication). 

In Fig. [6j we eventually analyse the dependence of the 

mean interface displacement (s) = \J ^ Sd^^ (^) 
solid fraction fsi for a fixed filter length L^iff = 0.75 hmesh^ 
where hmesh is the characteristic mesh size. An increasing 
branching for larger solid fraction is clearly seen in the op- 
timal designs, which results in a better support for the inter- 
face. In the chart we compare the mean displacement for the 
optimal design to the displacement for a post of circular cross 
section and same fsi. We can see that the optimised design 
always performs better than the simple circular cross section, 
and even more so for large solid fractions, which is again a 
consequence of the higher degree of branching in the opti- 
mised configuration. 

Conclusion and outlook 

In this paper we applied topology optimisation to the stability 
of superhydrophobic surfaces. We found that this technique is 
very effective for the task. Branching structures are found to 
be optimal to support hydrostatic pressure for a Cassie-Baxter 
state, in a two dimensional analogy to natural structures. We 
also analysed the effect of a solid fraction constraint on the op- 
timal design, as well as the use of a PDE filter to obtain designs 
suitable for fabrication. Further work will include the fabrica- 
tion and characterization of such optimised microtextured sur- 
faces. Preliminary fabrication results obtained at DTU Nan- 
otech suggest that the optimal shapes can be reproduced with 
a high degree of precision using common lithographic tech- 
niques. A further step will be to use a cost effective procedure, 
such as injection moulding, to produce the same designs. 

This research is funded by the NanoVation consortium. The 
authors thank Kristian E. Jensen and Rafael J. Taboryski for 
useful suggestions and discussions. 
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